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1. Introduction 



The fundamental quantities energy-momentum, and angular momentum of the 
gravitational field are elusive. Globally, for spacetimes which are asymptotically flat 
(or even anti-DeSitter), there are well defined values for the total energy, etc., given 
by the Bondi and ADM expressions integrated over spheres at null infinity and spatial 
infinity, values which are directly related to quantities that can be physically measured 
by distant observers. However, unlike all other fields, for the gravitational field these 
quantities have no well defined local density since, according to the equivalence principle, 
an observer cannot detect any features of the gravitational field at a point. It has even 
been argued that the proper energy-momentum of the gravitational field is only total; 
that it cannot (or should not) be localized (see, e.g., |J p 467). 

Localization is certainly possible if it is simply understood to mean "find some 
way of dividing up the total". In particular each of the many proposed pseudotensors 
provides for such a localization. But pseudotensors depend on the reference system, so 
they actually provide for many different localizations each of which includes some rather 
arbitrary unphysical content; thus one can arrange, at any selected point, for almost 
any value: positive, zero or even negative. Certain positive energy proofs (e.g., 0, f|) 
may be a better alternative: each gives localizations which are positive. But actually 
these "positive localizations" are not truly "local" , for they divide up the total in a way 
which depends on the configuration nonlocally, since they each depend on the solution 
of an elliptic equation which in turn depends on the field values everywhere. 

Appreciating the fundamental nonlocality of the gravitational interaction yet 
believing in the basically local nature of physical interactions led to the idea of quasilocal 
quantities: quantities that take on values associated with a compact orientable spatial 
2-surface. Expressions for quasilocal quantities in the context of Einstein's general 
relativity theory have been proposed from many perspectives including null rays ||, 
twistors ||, a fixed background 0, symplectic reduction ||, spinors || 10, 11], |i"2 , a 2+2 
"Hamiltonian" fll3 |, Hamilton- Jacobi |l4j and Ashtekar variables [0. Not surprisingly 
the various definitions generally give different results [16|, 10] . 

Lists of criteria to be satisfied by a quasilocal energy have been devised, see, e.g., 



17, OT. Usually it is required that it should vanish for flat spacetime, give reasonable 



values for weak fields and spherically symmetric solutions, and should approach the 
ADM and Bondi values in the appropriate limits. Opinions differ concerning whether 
quasilocal energy must be non-negative, but in any case such criteria are not sufficiently 
restrictive — in fact there remain an infinite number of possibilities ||T0 |. 

Here we offer a comprehensive presentation of some new ideas regarding quasilocal 
energy, many of which were first developed in |18[ and briefly reported in ||19|| . In our 
approach we begin with the idea that energy is naturally associated with time translation 
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and is thus given by the value of the time translation generator: the Hamiltonian. Hence 
good expressions for quasilocal quantities (energy, etc.) should be based not only on the 
variational principle and Noether's theorem but especially on the canonical Hamiltonian 



procedure, a viewpoint shared by some other investigators, e.g., [§, |14 

The fundamental feature that distinguishes our Hamiltonian based expressions for 
quasilocal quantities (aside from the fact that we treat rather general geometric gravity 
theories) is that they are 4-dimensionally covariant. We believe this is an essential 
quality for a physically meaningful (observer independent) definition of quasilocal 
quantities appropriate to a covariant gravity theory. Technically, our covariant 
Hamiltonian formulation p0[ is obtained by using differential form techniques J21 



in 



lieu of the usual spacetime splitting with its associated loss of manifest 4-covariance. 

For dynamic geometry theories the Hamiltonian has a special form connected with 
its role as the generator of displacements along a timelike vector field N. Noether's 
theorem applied to local translations along a vector field reveals that the Hamiltonian 
3-form (density) has the special form 7i(N) = (terms proportional to field equations) 
+dB(N). Consequently, the Hamiltonian H(N) = Jj,7i(N), which displaces a finite 
spacelike region E along N, has a value (on a solution to the field equations) given 
just by ;f> dT; B(N), the integral of the 2-form B(N) over the boundary of the spatial 
region. Thus this boundary integral will determine the value of the quasilocal quantities. 
Although the Hamiltonian field equations themselves do not depend on the boundary 
term, the expression for B(N) (unlike the case for other Noether conserved currents) is 
nevertheless restricted by the Hamiltonian variational principle. 

To fix the Hamiltonian boundary term we consider the variation of the Hamiltonian 
density. In general we get an expression of the form 57i(N) = (field equation terms) 
+dC(N). The total differential gives rise to a boundary integral term in the variation 
of the Hamiltonian. As Regge and Teitelboim have nicely explained p2f , it is necessary 
that the boundary term in the variation of the Hamiltonian vanish (only then are the 
functional derivatives well defined). This is not a problem for finite regions — if we 
fix the appropriate quantities on the boundary. However, when we consider the limit 
r — > oo,the boundary term § C(N) for gravity theories does not vanish asymptotically 
in general. (In particular this is in fact the case for Einstein's theory.) To compensate 
for this the B(N) term in the Hamiltonian needs to be adjusted. In this way the form 
of B(N) at infinity is constrained and the value of the total conserved quantities fixed. 

This argument was applied to the Poincare Gauge theory (PGT) and even more 
general gravitational theories for both asymptotically flat spaces and asymptotically 
constant curvature spaces by one of us [EJ], [23| with important improvements by Hecht 



. Expressions for the Hamiltonian boundary integrand which give the correct total 



conserved energy, momentum and angular momentum at spatial infinity |25] and at 
future null infinity [EHJ were obtained. These same boundary expressions can also 
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be applied to a finite region to give quasilocal values. But the very form of Hecht's 
improvement helped to show the way to other equally valid expressions. 

A comprehensive and systematic investigation of the role of total differential terms 
(i.e., boundary terms) in both the Lagrangian and Hamiltonian variational principle is 
really needed. Years ago Kijowski emphasized the importance of symplectic methods 
pTJ in such investigations. More recently there has been increasing recognition of the 
importance of symplectic methods and of boundary terms in variational principles; we 
have benefited from the progress made by several workers 00, 130, 13T], B3, 133], [T3 . 



A nice mathematical theory has been developed, which purports to incorporate the 
fundamental principle of physical interactions, and is especially well adapted to our 
problem. This is the theory of symplectic relations first proposed by Tulczyjew []34f 



and developed with Kijowski [j35| . (Later Kijowski and coworkers applied it to gravity 
theories, e.g. [§, |36fl . More recently Wald |37j has independently obtained some nice 
symplectic results.) In this treatment the (normally neglected) boundary term in the 
variational principle reflects the symplectic structure which, in turn reveals what physical 
variables are the "control" variables (held fixed) and which are the "response" variables 
(determined by the physical system). 

We shall use this theory to better understand both the general dynamic geometry 
Lagrangian variational principle and its associated boundary terms and the covariant 
Hamiltonian formulation. Hence the symplectic formulation is a key principle underlying 
our work. With its aid we will recognize the symplectic structure of the boundary term 
C(N) in the variation of the Hamiltonian and be guided to select those Hamiltonian 
boundary terms B(N) which give rise to a "covariant" symplectic structure for C(N). 

The application of these ideas, the covariant Hamiltonian procedure along with a 
covariant symplectic boundary variation, determines the expressions for the quasilocal 
quantities of a gravitational system. In this way we obtain a Hamiltonian 3-form 7i(N), 
which generates the evolution along the vector field N, and which includes, for each 
independent variable (e.g., the connection, the coframe, the metric), one of two possible 
covariant boundary terms depending whether the field or its conjugate momenta is held 
fixed (controlled) on the boundary. Then the variation of the Hamiltonian, in addition to 
the field equations, includes a boundary term with a covariant symplectic structure which 
reflects the choice of control mode. For each geometric field (metric, frame, connection) 
we can control either the field value or the momentum (essentially, Dirichlet or Neumann 
boundary conditions). Each choice gives rise to a different quasilocal value. Thus, as in 
thermodynamics, there are several different kinds of "energy", each corresponds to the 
work done in a different (ideal) physical process. 

The Hamiltonian boundary terms are the expressions for the quasilocal quantities. 
They determine the value of the Hamiltonian. The physical meaning depends on the 
displacement: energy for a time translation, momentum for a space translation, and 
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angular momentum for a rotation. (Although our Hamiltonian formalism apparently 
presumes a timelike displacement, this is really no limitation, for the Hamiltonian 
is linear in the displacement vector field. Hence considering the difference between 
two suitably chosen timelike displacements gives the meaning of the Hamiltonian for 
spacelike displacements such as translations and rotations.) 

A noteworthy feature of our expressions is that they have an explicit dependence 
on a selected reference configuration. The reference configuration has a simple meaning: 
when the dynamic fields have the reference configuration values on the boundary then 
all quasilocal quantities vanish. In particular, the reference configuration determines 
the "zero" of energy. 

We present the formalism for rather general geometric gravity theories. Specifically 
gravity theories which can have an independent metric and connection. The connection 
need not be metric compatible nor symmetric. The field equations are presumed to 
follow from a Lagrangian which may depend in any way on the metric, the curvature, the 
torsion and non-metricity tensors (but not on their derivatives). This general geometric 
approach has proved to be a good guide for specific and more specialized theories. Our 
expressions are easily restricted to special cases such as Riemann-Cartan, Riemannian 
or teleparallel geometries. For Einstein's theory and the Poincare gauge theory (for 
asymptotically flat or asymptotically constant curvature solutions) our expressions are 
related to previously known expressions (e.g., [^2[ [39|, |40], [41] ) and give good values 



for the total conserved quantities at spatial |25|] and future null infinity [^6[ . 

The outline of this work is: in section 2, we give the covariant Hamiltonian analysis of 
a general geometric gravitational theory in the language of differential forms. In section 
3 our requirements and expressions, based on the covariant-symplectic Hamiltonian, for 
quasilocal quantities are presented. We apply the ideas to black hole thermodynamics 
obtaining the generalized first law and an expression for entropy for these general 
theories in section 4. The specialization of our expressions to certain gravity theories 
is discussed in section 5. We restrict our relations to general relativity, apply them to 
spherically symmetric solutions and compare them with some other authors' results in 
section 6. Finally our concluding discussion considers the various choices in selecting a 
Hamiltonian boundary term-quasilocal expression, noting in particular how this scheme 
provides for an orderly system relating the choices to physically meaningful properties 
on the boundary of a spatial region. 



2. Covariant Hamiltonian formalism 

In this section we briefly recount the relevant features of the covariant Hamiltonian 
formalism |2(], |23| . We consider quite general dynamic geometry theories of the metric- 
affine type (see, e.g., jEJ) as well as some important special cases. The possible 
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geometric potentials are the metric coefficients g^ v the coframe one-form $ a and the 
connection one-form u a p. The corresponding field strengths are the non-metricity one- 
form 

Dg^ := dg^ - uj 1 ^g lv - uP v g^, (1) 
the torsion 2-form 

0" -=Dd a :=d$ a + oj a p A$P, (2) 
and the curvature 2-form 

n a fs := du a p + A ufp. (3) 

The usual approach is that (second order) dynamical equations are presumed to 
follow from a variational principle. We presume that the Lagrangian density is a 
scalar valued 4-form depending only on the gauge potentials and their first differentials: 
£(2) _ C^ 2 \g,'&,u;dg,d'&,du}). (Note that we have, in a very natural way, excluded 
higher derivatives of the fields. However, there is really no obstacle in principle to 
extending our formalism to include them, say via the expense of introducing extra fields 
and Lagrange multipliers.) In view of the invariance under local Lorentz transformations 
the more covariant form: C^ 2 \g, Dg, O, Q) has advantages. Instead of this usual 2nd 
order type Lagrangian 4-form we use a "first order" Lagrangian 4-form: 

C := Dg^An^ + Q a Ar a + Q a pApJ - A(g, vr, r, p). (4) 

Independent variation with respect to the potentials (g,<d,uj) and their conjugate 
momenta (tt, r, p) yields 

.„ . 5C cnn 5C . „ 5C 5C . .... 5C 

+ r^a A + ^l**"-""" + W° A T„ + fo^ A /)„"), (5) 
op a 

which implicitly defines the first order equations (their explicit form will not be needed). 

The Hamiltonian formulation displays the "time" evolution of the physical system. 
To proceed in a covariant way we consider instants of "time" to be given by spacelike 
hypersurfaces determined by a function t. A "time evolution" vector field N displaces 
one constant t space- like hypersurface S t to another; consequently i^dt = 1. Time 
evolution is then given by the Lie derivative where £n = iNd + diw on the components 
of form fields. We regard t and N as being fixed when we vary our physical variables. 
One way to obtain a set of Hamiltonian equations is the "spatial restriction" (pullback) 
and "time projection" (evaluation on the vector field N) of the first order equations. This 
would give the initial value constraint and dynamical evolution equations respectively. 
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The same equations can also be obtained from a Hamiltonian. To obtain the 
Hamiltonian in a covariant form we decompose the Lagrangian as follows: 

C = dt A ijqC 

= dtA (£ N 9^ U + £ N $ a Ar a + £ N u a p A pj - H(N)). (6) 

In this way we find the Hamiltonian 3-form 7i(N) = N^Ti^ + dB(N), where 

N»H„ := i N A + Dg^ A % N ^ V - Q a A i N r Q - Q a p A i N pJ - i N $ a A Dr a 

- i N u a p (DpJ - g a y v - g m ^ +0t>A r a ), (7) 
B(N) :=i N d a T a + i N u a p P J. (8) 

The Hamiltonian 3-form 7i{N) is the same quantity which arises when one considers 
Noether's second theorem for a local translation along a vector field N applied to the 
first order Lagrangian (f|), as we now briefly explain. 

Geometric theories are invariant under local diffeomorphisms. For a first order 
Lagrangian of the form £ = dip Ap — A the general variational formula 

5C = d(5cp A p) + 6<p A — + — A 6p, (9) 

dip bp 

must be identically satisfied if the variations are generated by an infinitesimal 
diffeomorphism. Hence 

£nC = dixC- = d(£]y{p Ap) + £n<P A h — A £nP- (10) 

dip op 

Consequently the "Hamiltonian 3-form" 

H{N) := £ N p Ap-i N C, (11) 
satisfies a differential identity 

dH(N) = (terms proportional to field equations). (12) 

In addition to this Noether differential identity, because displacement along iV is a 
local symmetry, we also have an algebraic identity. Using £ = dp A p — A in ([□]) we 
rearrange the Hamiltonian into the form 7i(N) = N^Ti^ + dB(N) and then substitute 
dH(N) = d{N fX H ll ) = dN^ A + N>*dH» into the differential identity flTJ). The 
coefficient of dN^ gives an algebraic identity, 

7Y M = (terms proportional to field equations). (13) 

From (0) we see that "on shell" (i.e., when the field equations are satisfied) the 
Hamiltonian 3-form is a "conserved current" and the value of the Hamiltonian H(N), 
the integral of the H(N) over a space-like hypersurface E, is a conserved quantity. 
From the algebraic identity (|T3"D, vanishes when evaluated for solutions of the field 
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equations. Thus, "on shell", the value of Hamiltonian H(N) depends only on the term 
dB(N). By the generalized Stokes theorem, the value of the Hamiltonian is just the 
integral of the 2-form B(N) over the 2-surface boundary <9E of the 3-hyper surf ace E. 
This value should determine the energy and the other quasilocal quantities. However, 
just as with other Noether currents, TC(N) is not unique; we can add to it a total 
differential without affecting the conservation relation fll2|), this amounts to modifying 
B(N) (8). Such a modification does not affect the Hamiltonian equations of motion. 
The Hamiltonian field equations are obtained by varying the action (6). The variation 
of the above Hamiltonian has the form 

5H(N) = (field equation terms) + dC(N), (14) 

where 

C(N) := i N (5g^ + 5$ a Ar a + 5uj a p A pj). (15) 

A modification of B(N) (8) affects only the total derivative term dC(N) which, upon 
integration over E becomes a boundary term. 



3. Expressions for quasilocal quantities 

Our purpose is to find quasilocal quantities for gravity theories. Hence we want to 
determine the expression for the term B(N) in the Hamiltonian. From the variational 
principle we know that a total differential term can be added to the Hamiltonian without 
changing the field equations. But this, of course, will change the expression of B(N), and 
thereby the value of the quasilocal quantities. What is the proper "physical" expression 
for quasilocal quantities? We have already noted that the requirement of having the 
desired limiting values is far too weak. We use certain additional theoretical criteria 
involving the boundary term in the variation of the Hamiltonian to severely restrict the 
form of a "good" expression for the term B(N). 



3.1. Requirements 

Consider the total differential term C in the variation of the Hamiltonian. We propose 
three theoretical requirements on the term C in order to limit the quasilocal boundary 
expression B: 

• [Rl] Well-defined Requirement 

— The term C must vanish at the boundary in general both for finite and infinite 
regions. 

• [R2] Symplectic Structure Requirement 

— The term C must reflect the symplectic structure $5], which tells us which 
physical variables are "control" variables (i.e., we control them at the 
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boundary) , and which are "response" variables (i.e., are determined by physical 
laws). 

• [R-3] Covariant Requirement 

— The control and response variables must each appear in the form of a covariant 
combination (as befits a covariant theory). 

With respect to the Well-defined Requirement Rl, the expression of the term C(N) in 



equation (JL5|) is no problem for a finite boundary. We need merely fix the appropriate 
variables to vanish at the boundary. For a boundary at infinity, however, one needs 
to consider a limit and the natural rate of fall off in the asymptotic region. In this 
case, the term C(N) in equation (|15D is nonvanishing generally. (Specific examples are 
the Einstein theory |22|] and the Poincare gauge theory (PGT) for solutions which have 
asymptotically zero or constant curvature.) In order to correct C(N) to satisfy Rl, one 
must modify the term B in equation (8). 

3.2. Well-defined version 

Some time ago within the context of metric compatible theories, i.e., the PGT, one of 
us |20| proposed a modified version of the boundary term (8): 

B(N) = i N $ a Ar a + i N oj a p ApJ 3 + Au a p A i N pJ. (16) 

With this adjustment the total differential term C(N) in the variation of 7i(N) takes 
the form 

C(N) = 5i N $ a AT a - 5t) a A i N r a + 8i N uj a p ApJ + Au a p A 8i N pJ, (17) 

where, for any quantity a, a is the reference configuration value of a which is fixed 
under variation, and Aa := a — a. This expression should vanish at spatial infinity, in 
particular for asymptotically flat or constant curvature PGT solutions. But for some 
theories, e.g., the PGT with a massless torsion field P3"f , the boundary term B(N) 
still needs adjustment. Hecht recognized this and proposed an important modification 



24 1 to the i^uj factor. (We have already noted that the Hecht expression has since 
been successfully tested at spatial and null infinity for both asymptotically flat and 
asymptotically constant (negative) curvature solutions p6[|.) 

However we can also add a term Ai^^a into the original expression fllTf) to give 
the formula more symmetry. Then 

B'(N) = i N i) a AT a + At) a A i N r a + i N oj a p ApJ + A^p A i NP J, (18) 

and 

C\N) = 5i N $ a Ar a + Ar) a A 5i N r a + 5i N u a p ApJ + Acu a p A 5i N pJ. (19) 
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From these considerations we see that the Weil-Defined requirement Rl cannot 
determine the quasilocal term B(N) uniquely. Moreover neither C(N) nor C'(N) satisfy 
the Covariant Requirement R3. 



3.3. Covariant symplectic structure version 

More recently we have improved these expressions and obtained several different 
boundary terms, all of them satisfy the requirements Rl, R2 and R3. They include, for 
each independent k-form field <p (e.g., the connection, the coframe, the metric), one of 
two possible covariant boundary terms: 

B V (N) = i N ipAAp-(-l) k AipAi N p or B P (N) = i N (pAAp — (—l) k AipAi N p, (20) 

depending upon whether (p or its conjugate momenta p is held fixed (controlled) on the 
boundary. Then the variation of the Hamiltonian, in addition to the field equations, 
includes a boundary term which is a projection onto the boundary of a covariant 
symplectic structure: 

di N (5(p A Ap) or di N (— A(p A Sp), (21) 

respectively, which reflects the choice of control mode. 
Thus for the geometric fields we take 



B(N) 



-Ag^ u i N > v 1 f %N^Ar a + A$ a A i N r a 
-Ag^ u i N n^ J 1 i N i a Ar a + A$ a A i N r a 



- \ , o„ a „ ft , A„,a a • rn {5 ' 



i N u a pApJ + Au a pAi N ° Pa P 
i N u a i3 Apa 13 + Au a f3 A i Npc 

where the upper (lower) line in each bracket is to be selected if the field (momentum) 
is controlled. Hence, as in thermodynamics, there are several kinds of "energy", each 
corresponds to the work done in a different (ideal) physical process [pSJ, 0. 

For the geometric variables, the total differential term in 87i{N) is of the boundary 
projection form di^C with C now given by the covariant expression 

J 5g, u An^ If 5"d a A Ar a I f 5u a A Ap/ 1 

\ -Ag^ u Stf" J + 1 - A 5r a J \ -Aa;^ A 5pJ J ' 1 j 

where again the upper (lower) line in each bracket corresponds to controlling the field 
(momentum). 

Of course we cannot expect to get completely covariant expressions when we 
have a non-covariant dynamic variable such as a connection. Note however that Au, 
being the difference between two connections, is tensorial, so our quasilocal boundary 
expressions are covariant — aside from the manifestly non-covariant explicit connection 
terms in B. These i^u terms include a real covariant physical effect plus an unphysical 
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(noncovariant) dynamical reference frame effect. These two effects can be separated by 
using the identity 

(i N u a p)^ = i N Q a + DN a - £ N $ a , (24) 

to replace the znuj factor within the last bracket in (22) with two covariant terms plus 
a manifestly non-covariant term of the form (i^t? ) pApJ 3 '. This £/v$Ap piece is really 
needed in our general Hamiltonian to generate the dynamics due to the gauge freedom of 
the reference frame — but it represents an unphysical (observer dependent) contribution 
to the quasilocal quantities. For the purposes of calculating the physical value of the 
quasilocal quantities it should be dropped. (Note that if N is a Killing vector then this 
term can be made to vanish for a suitable choice of frame.) 

3.4- Uniqueness 

An important property of our expressions is uniqueness. Of course, the uniqueness 
property depends on the proposed requirements. Under our three requirements (Rl, R2 
and R3), our expressions are the only possible ones. 

According to our covariant symplectic structure requirement (R2 + R3), for each 
field potential tp and its momenta p, the boundary term of the variation of the 
Hamiltonian must be one of the following two types 

di N (5(pAAp) or di N (-Aip A 8p). (25) 

The only acceptable modification to the Hamiltonian which preserves the covariant 
symplectic structure requirement is of the form di^J- *, a projection of a 4-covariant 
form onto the boundary. The 3-form T must depend on the control variables only 
algebraically f and must be scalar-valued. Working with our dynamic variable pairs 
(sV' 7r '"0> {® a i T a) an d (^ Q /3, PoP), the only acceptable combination is of the form 
di N (A(p A Ap). This just switches between the two control modes: 

B p = B v -i N (A(pAAp). (26) 

Of course one could add some covariant quantity which doesn't depend on the dynamic 
variable. That would leave the Hamiltonian variation symplectic structure unchanged 
but would "renormalize" the quasilocal quantities — albeit in an "unphysical way". We 
rule this out by normalizing our quasilocal expressions so that they vanish if the dynamic 
variables equal the reference values on the boundary. 

Although we have determined the expressions 'uniquely', there still exist two 
undetermined things. One is the displacement vector field N. For example, if we 

f Since if it depends on the differential of the control variable, we must fix the differential of this 
variable at the boundary. This is forbidden by the symplectic structure. 
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want to calculate the energy of a physical system in a finite region, we have not yet 
specified how to select the appropriate timelike vector N. The other is the reference 
configuration. We will look closer at these issues later. 

3.5. Spacetime version 

Now let's go back to the Lagrangian level using the same requirements in order to get a 
completely covariant analysis. In equation (5), the boundary term satisfies the Covariant 
Symplectic Structure Requirement R2 and R3, but not the Well-defined Requirement 
Rl. 

For the purpose of getting boundary terms satisfying our requirements, we 
considered ways to modify the Lagrangian by adding an extra total differential term. 
In a first order Lagrangian we found that we could include a boundary term by letting 

£ := dtp A p - A - dlC, (27) 

where K = Atp A p or /C = Aip A p. Then 

5C = (field equations) + d( j ( 28 ) 

We can relate this to the Hamiltonian analysis for 

5in£ = %nSC = ^(field equations) +ijyd( j ^ ^ f jV (29) 



-Aip A Sp 



The latter term can be rearranged into 



which consists of a total time derivative (which integrates into a term at the initial and 
final times and constitutes a canonical transformation) and a spatial boundary term 
which is identical to our covariant symplectic Hamiltonian boundary variation term. 
By the way, it is here that one can see the covariant origin of our required form for the 
projection of a covariant quantity fl2lf) . 



4. Black hole thermodynamics 



Here we will use arguments similar to those of Brown and York [44, 45 1 and of Wald 



see also ||7J ) along with our symplectic Hamiltonian expressions to derive the first 
law of black hole thermodynamics for these general geometric gravity theories. 

For the theories of gravity that we are considering, the null rays and causal structure 
are governed only by the metric and its associated Riemannian geometry. Hence the 
basic geometry of black holes in these theories should be the same as that of GR. 



13 



The quantity corresponding to the temperature in black hole thermodynamics is the 



surface gravity, k, defined at any point of a Killing horizon by [48 

fV ft /=: KX ' 3 1 (31) 

where \ a is the Killing vector field normal to the Killing horizon and V is the 
Riemannian covariant derivative. Hawking proved that the event horizon of a stationary 
black hole is a Killing horizon in general relativity fl49 |. For certain general gravity 



theories it has been shown that if a black hole has a bifurcate Killing horizon (the 
Killing vector vanishes at this horizon) then the constancy of the surface gravity holds 
||50|| . For the theories under consideration here we presume that the black hole has a 
bifurcate Killing horizon and that the zeroth law is satisfied: the surface gravity (which 
is proportional to the temperature) is constant over the entire horizon [|51[ . 

The first law of thermodynamics concerns the conservation of energy. Because 
our Hamiltonians satisfy a conservation law, the corresponding formula for gravitating 
systems, in particular black holes, can be obtained and the quantity corresponding to 
the entropy identified. 

The basic idea is to consider the region between the horizon and infinity. We wish 
to vary certain parameters such as the total energy and angular momentum. However in 
the formulations discussed earlier such quantities were not control parameters. To obtain 
the necessary additional free variations on the boundary we must allow the displacement 
vector field to vary. Then the variation of the Hamiltonian 3-form, including the 
variation of the displacement vector field, is of the form 

5H(N) = (field equation terms) + (5N) fl n f , 

+ <i(iAr(symplectic control mode terms) + B(SN)) , (32) 

with the 7f M term vanishing by the initial value constraints. Note that the response 
variable for the displacement vector is the "component" of the quasilocal quantities. 
Here B can be any one of our boundary terms representing any control mode as long as 
the "reference frame energy" is removed by dropping the £n$ contribution via the i^uj 



replacement (|24f) . Following Brown and York [f44] , we now introduce the microcanonical 
Hamiltonian via a Legendre transformation: 

H^ CXQ {N)=H{N)-dB{N). (33) 

The variation of the microcanonical action is of the form 

<KSmicro[iV]) = J (field equation terms) - (SN^H^ 

— f (ijv(symplectic control mode terms) + (5B)(N)) . (34) 

The variation of .Smicro should vanish on a solution with the control variables fixed at 
the boundary. Suppose the hypersurface includes a boundary with two components, 
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one at infinity and the other at the bifurcate Killing horizon, H. Then we get a kind of 
energy-momentum conservation law: 

= I (5B)(N) = £ {5B){N) - I (5B)(N) . (35) 

JdT: Joo Jh 

Now we shall choose our variations to be such that they perturb the asymptotic values 
of the total quantities. 

Let X t and X v denote the Killing vector fields which approach the time translation 
and rotation at the infinity. Then the total energy is the value of H(X t ) = ^ B(X t ) and 
the total angular momentum is the value of H(X^) = ^ B{X tp ). It is worth remarking 
that this is true for all the possible control modes, asymptotically they give the same 
total results for they differ by terms with more rapid fall off. 

We assume as usual that the Killing vector x can be decomposed into 

X = X t + £l H X v , (36) 

and vanishes on the bifurcate Killing horizon, where Qh is the "angular velocity of the 
horizon" . 

Now let the displacement N be the Killing vector x- Then 

/ {SB) (N) =5E + n H 5 J. (37) 

J oo 

On the other hand, for the surface integration over the horizon boundary, because the 

o 

Killing field vanishes on the horizon the only terms that can contribute involve Dx or 

o 

Dx = Dx + Ac^x, both of which which reduce to Vx on the horizon (since the affine 
connections differ from the Levi-Civita connection by tensor terms multiplied by the 
vector field which vanishes on H). Hence we get 

/ (5B)(x) = I Vpx a $P« P = k I e p a 5pj =: kSS, (38) 

here we have assumed that V a Xp = Ke a/3, where e a p is the bi-normal to the horizon. This 
has the form of the first law and identifies an expression for the black hole entropy S for 
our general gravity theories to be dependent on the horizon components of p^, which 
is the field momentum 2-form conjugate to the curvature. This result is consistent with 
the well known Einstein theory value, where the integral of e a ^p a/ g is proportional to 
the area. Remarkably, it is of the same form as an expression found previously for other 
general gravity theories which were based on Riemannian geometry with no torsion 
or non-metricity but, on the other hand, were allowed to contained arbitrarily high 
derivatives of the curvature |5T|, f25|. These works used various techniques ranging 
from the Noether- Charge analysis of Iyer and Wald to the microcanonical functional 
integral method of Brown. The extent to which these results can be combined with our 
work is not yet clear to us. 
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5. Some theories of gravity 

Our general formalism readily specializes to either coordinate or orthonormal frames, 
moreover it is suitable for the most general metric-affine gravity theories H2~, 



(For earlier proposals for quasilocal quantities for such theories see §5.7 in 
Nevertheless, it also readily specializes to less general geometries like Riemann-Cartan, 
Riemannian or teleparallel and to gravity theories formulated in such geometries. 



Teleparallel theories (i.e., curvature vanishes), see, e.g., |53|, [54j], are somewhat 
special; hence we leave the treatment of their covariant Hamiltonian formalism, and 
quasilocal quantities for future work. For the Poincare gauge theory (PGT) J55], 
one could impose the metric compatible connection condition via a Lagrange multiplier. 
However, a neater method is to use orthonormal frames and restrict the connection 
algebraically: u al3 = —u 13 ". Then one need merely drop all the g, n terms from our 
expressions. The most general potential A($, r, p) contains 3 independent quadratic 
terms in r and 6 independent quadratic terms in p: 

r a = --*(E«n (n) 6 Q ), p^ = -^-i*(^&n (n) ^), (39) 

which express the momenta linearly in terms of the algebraically irreducible parts 
of the torsion and curvature. Here we have introduced the convenient notation 
rf- : = *($ a A • ■ •), with 77 := *1 = $° A tf 1 A $ 2 A $ 3 . 

The Einstein-Cartan theory is a special case of the PGT with the Lagrangian 

2% X 

(where A co <, is the cosmological constant and % := 8tcG/c a = 8tc in geometric units). 
The variables are the orthonormal frame and the metric compatible (antisymmetric) 
connection one-form. The associated conjugate momenta are 

r a = 0, p J = L ria P. (41) 
2% 



-EC = i- + (40) 



Properly we should introduce these momenta and regard ([|l]) as constraints which can 
be enforced with Lagrange multipliers. However, because these "momenta constraints" 
are purely algebraic, we can take a short cut and directly use (ffl) as our first order 
Lagrangian. The corresponding quasilocal expressions are 

B u = ^-(Auj a fs A i N riJ + %n^p At?/) , (42) 
^X 

B# = ^-(Au a p A i NV / + i N u a p A77/) . (43) 
^X 

The gravity theory which attracts the most interest is Einstein's General Relativity. 
The Riemannian geometry of GR does not have an independent connection while our 
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formalism is geared to independent variations of the frame, metric and connection, Yet 
there is no incompatibility, indeed one can proceed in many ways including: (1) impose 
the metric compatible and vanishing torsion conditions via Lagrange multipliers, (2) 
use coordinate frames, algebraically impose the symmetry of the connection coefficients 
and impose metric compatibility via a Lagrange multiplier, (3) use coordinate frames, 
algebraically impose the symmetry of the connection coefficients and (for vacuum and 
all non-derivative coupled sources) obtain the metric compatibility condition as a field 
equation via a Palatini type variation, (4) use the orthonormal frame EC theory 
and impose the vanishing torsion condition via a Lagrange multiplier, (5) use the 
orthonormal frame EC theory and obtain (for vacuum and all non-derivative coupled 
sources) the vanishing torsion condition as a field equation. 

The latter method is technically the simplest and we will use it here to obtain 
our quasilocal quantities for GR. All of the other procedures give similar results. 
Similar but not necessarily identical — because the boundary control differs — for example 
Aw"^ A 5r] a P = Au a /3 A <5i9 M A rjj 5 ^ in terms of orthonormal (co)frames, but it 
equals Ar^ 5{^gg p<7 ) 8 v a l d 3 x„ in a coordinate basis. A comprehensive investigation 
covering all these options is underway. 



6. Quasilocal values in general relativity 

We now consider our quasilocal expressions for general relativity. This will help to give 
a better understanding of their mathematical nature and physical meaning and will 
permit a comparison with the results obtained by other authors. 

The Hilbert Lagrangian fl40|) of the Einstein(-Cartan) theory (dropping the 
cosmological constant and often the constant 2x for simplicity) can be spacetime 
decomposed according to 

£ = dt A zat£ 

= dt A [£ N u% A rjj + N»n a p A V J '„ + i N uo a p D V J - d(i N u a p V J)] , (44) 

which identifies the Hamiltonian 3-form as 

H = -N»n a p A T] a \ - i N u a p D V J + d{i N u a p V J) . (45) 

On a constant t = x° surface, the principal term in this Hamiltonian, —N^Ci, 01 p Ar) a P ^ = 
N^2G U fj^v, just reduces to the familiar ADM value 2N^G° fl (the rotation gauge 
generator term is proportional to the torsion). However the boundary term still needs 
adjustment [22]. We propose using one of the Hamiltonian boundary terms (42, 43). 
Removing the observers "reference frame energy-momentum" ~ £n$At) via the identity 
([24]) gives our covariant quasilocal boundary expressions for GR: 

£L = Au a p A i N ?iJ + DpN a At?/ , (46) 
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B# = Auj a p A i N r}J + D p N a ArjJ 



(47) 



For the simple alternative choice of a coordinate basis with the metric and 
symmetric connection as the variables, the derivation of the Hamiltonian and its 
boundary terms is formally the same. The control variables are then ^/—gg^ a and 
T a p\ respectively, and the components of the quasilocal expressions (46, 47) take (with 
B := {l/A)B T Pe Tpixv dx^ A dx u ) the form 



BY = N p (-g)^ <f Ar> + DpN a A(^g^) 5% , 

,8a \ rrp 



9 

B T / = N»(-g)-2 g^ Ar a p\ + D p N a A 



-gg 



(48) 
(49) 



These new quasilocal boundary expressions are similar to certain results obtained 
by earlier investigators. With a unit displacement, the first term of B g by itself, 



U^(g, AT) := (-0)3 gP° AT a 0x ^ > 
if we use a vanishing reference connection, reduces to the Freud superpotential |5 



U, rp (g,T) 



-g) 2 g^d x \ 



-g(g TU g pX 



g pv g rX )} ■ 



(50) 



(51) 



Using just this superpotential as the Hamiltonian boundary term is equivalent to 
integrating its divergence, the Einstein pseudotensor, over the spacelike hypersurface. 
Our second term, without the A, has the form 2D^ P N T \ which is the well known Komar 



expression |pq| . Several investigators, working from different perspectives, have arrived 
at an expression which is the sum of these two terms: 



N^(g,AT) 



^f=gg^ T 5 p] b„N x = 2\f—Q N^K^ + 2 
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-g 



D [p N 7 



(52) 



where K p := 2g /3 ' A Ar p ' p\. This was obtained by Chrusciel [59] and also by Katz 
using a background metric while Sorkin |60| used only a background connection. 
Chrusciel' s analysis is especially relevant to us since he also used the Kijowski-Tulczyjew 
symplectic-Hamiltonian ideas in much the same spirit as we have done; his work even 
includes a second expression corresponding to the conjugate control mode. Most of these 
earlier investigators recognized the quasilocal possibilities of their results, but they were 
primarily concerned with obtaining and discussing the virtues of this expression for total 
energy-momentum and angular momentum along with reconciling the Komar expression 
with those obtained by other techniques. Mention should also be made of an interesting 
covariant investigation of Ferraris and Francaviglia [31] which uses a global background 
symmetric connection and has also arrived at the same expression. 

Although the expression (^) obtained by the aforementioned investigators is quite 
close to our (48), an important difference lies in the presence of the A in the "Komar" 
term. In most practical quasilocal calculations this amounts to a 'difference which 
makes no difference'. However we believe that this A is really needed for the most 
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general considerations (in particular it is necessary to get the asymptotically correct 
behavior for the boundary term in the Hamiltonian variation). There are precedents for 



this factor beginning with the work of Beig and O Murchadha |38 



To see this it is necessary to do a space time split of our expressions. We '3+1' 
decompose expression (47) using an orthonormal frame with one leg e± normal to the 
hypersurface. The extrinsic curvature of the spatial hypersurface is related to the spatial 
restriction of certain 4-dimensional connection one-form components: 

-Kfltf , (53) 



while the spatial restriction of the spatial components u l j of the 4-dimensional 
connection one-form is just the 3-dimensional connection one-form (for details of the 
technique see |S2|). In this way we obtain 



B$ = Aoo ij A N x THjx + 2Aw ±j A N k rj ±jk + (D j N L - D^N j ) Ar} ±j , (54) 



where we have dropped certain terms, including DW^ Arjij, since they are proportional 
to I? -1 = N dt and thus vanish when restricted to the spatial hypersurface. Reverting to 
tensor notation with a coordinate basis, the first term is 

A^ a Vi j ± = ArVt^i = (g ij g kl - g ll g k] ) v l A 9k3 ^dS l , (55) 

where we have used 

AP W = y m (W k Ag ml + VAg mk - V m Ag kl ) . (56) 

Similarly, the second term can be rewritten in terms of the extrinsic curvature and then 
in terms of the ADM conjugate momentum: 

2Au^A N k V±jk = -2A(Ki t dx l ) A N k r] jk± = -2N k AK^ m Sfty/gdSi 

= -2N k A{K\ - Si&Jy/gdSi = -2JV fc A(-(/- 1 /W(/ J - fc )^d5i • (57) 
To appreciate the significance of the remaining term, we consider our boundary to be 

o 

at spatial infinity. Asymptotically, with constant in time and the usual fall-offs, the 
final term takes the form 

b°N x A$ k A Vxjk = g jl N tl \ g kn Ag m 5f k l Jg dS m ; (58) 

here we have made the acceptable assumption that (Ai) k ) rn is symmetric and hence 
determined by (1/2) Ag^. Taking it all together, at spatial infinity our boundary term 
reduces to 



i dSi 



Vg(g t3 g kl - g u g kj ) {NVAg kj + N ti A 9kj ) + 2Nm kl g . 



lk 



(59) 



where gij and V» are flat space quantities. This is just the Hamiltonian boundary term 
(generalized to a form suitable for computation in non-cartesian coordinates) found by 
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Beig and O Murchadha [38 to be necessary to allow for Poincare displacements of the 



asymptotic Minkowski structure at spatial infinity for the asymptotically flat Einstein 
theory. The novel feature which they discovered is the N^Agj^ term, which is, according 
to their analysis, absolutely necessary to account for the boosts. This is where the A 
in our version of the "M0ller-Komar" term affects the calculated values. Concerning 
the value of these boundary terms, the N and N k terms give the ADM energy and 
momentum respectively, the N k term also gives the angular momentum. Note that the 
N t k term makes no contribution unless the lapse N is asymptotically unbounded. 

More recently Katz, working with Bicak, Lerer and Lynden-Bell, has generalized 
his previously mentioned work, using Noether conservation arguments applied to a 
Lagrangian quadratic in the first derivatives of the metric along with a global (possibly 
curved) background, to obtain a boundary expression for conserved quantities which 
can be put in the very succinct form 

2{y/^N [p K T] + A{^^D [p N r] )}. (60) 



Via manipulations such as those used in conjuction with eq ( |52"D this is exactly the 
same as our expression B T g p . We regard this is a significant support for our particular 
GR results and thus indirectly for our general ideas (conversely, our independent 
approach also supports the work of these other investigators). It seems that our nearly 
simultaneous efforts are rather complimentary. Consequently the successes achieved by 
these other workers in applications as diverse as Mach's principle conservation laws 



at null infinity |53] and integral constraints for large cosmological perturbations are 
further conformations for our work. This interesting body of work has only recently 
come to our attention; unfortunately, perusal of the connections between our work and 
that of Katz et al and the applications suggested thereby cannot be developed further 
at this time. 

The work just mentioned, however, brings up an alternate approach to our whole 
analysis of quasilocal quantities. Instead of adjusting the Hamiltonian, we could have 
begun with an adjustment to the Hilbert Lagrangian by a boundary term. Our technique 
gives essentially only two different covariant versions for the GR Lagrangian with 
boundary term: 

C := n a p A - d(Au a p A riJ) , (61) 
C$ := Q a p A T]J - d(Au a p A r)J) . (62) 

Let expressions with a tilde denote the quantity with a trivial reference configuration, 
then 

Co := Wp A - d(u a p A rjj) . (63) 



The boundary term here is just the trace of the extrinsic curvature — this can be seen 
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by noting that the restriction of the boundary term to a surface with normal e± is 

2u ± a A V± a = -2K ab $ b A V± a = -2K ab g ab V± . (64) 

Hence we can identify as the Lagrangian used by York and Brown p9| (L4j]; its 
space-time decomposition has the form 

S 1 ■= [ C BY = ^- f <?x^g~R +- f d 3 xVhK - - [ d 3 xy/^Q , (65) 



hi 2xJm X J t> x J:i B 

where K and G are the traces of the extrinsic curvatures of the constant t and spatial 
boundary 3-surfaces which have metrics hij and 7^, respectively. Brown and York 
begin with this Lagrangian and, without discarding any total derivatives, construct 
the Hamiltonian along with its boundary term. Consequently, beginning with (^), a 
straightforward construction of the 'covariant Hamiltonian', without discarding any total 
derivatives, leads to their Hamiltonian boundary term expressed in terms of covariant 
quantities: 

B BY (N)=u a fi Ai N r } J. (66) 

By the same type of techniques as used above in connection with (|54"1,57), it is easily 
verified that this 4-covariant expression (it's essentially just the Freud superpotential 
again) really does yield the densities for the quasilocal energy and momentum fl4| which 
they derived from the physical Lagrangian S 1 . 

Brown and York took a Hamilton- Jacobi approach; they chose a reference action 
and considered S = S 1 — S°, so their renormalized quasilocal quantities are of the form 

B BY (N) := B BY (N) - B° BY (N) . (67) 

With a suitable choice for the reference action (it's necessary to use a reference 2- 
geometry which is isometric with the boundary 2-surface) we can obtain a covariant 
version of their expression: 

B BY (N) = u a p A i N r]J - u a fs A invJ = Acu a p A i N r]J. (68) 

Restricting to the 2-boundary, and using the previously referred to '3+1' decomposition 
techniques such as in (57) along with 

Auj^ A N\ j ± = 2N ± Aoo h A A r] h A ± = 2N ± A(-k AB $ B ) A r] h A ± 

= - 2N ± A(k AB a AB )r ]h± , (69) 

where ei_, <Jab are, respectively, the normal and metric (in an adapted frame) of the 
2-boundary, plus a '2+1' decomposition of (57), leads to the quasilocal energy and 
momentum densities of Brown and York: 

e = &lo , jt ■= -2(<7 tt (eH),II*7V/0 IS 1 , (70) 
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where k is the trace of the extrinsic curvature of the boundary 2-surface embedded in 
the constant t surface E. This covariant version of their expression is closely related to 
one of ours: 

B BY (N) = B#(N) - i N u a p AnJ . (71) 

They exactly agree only when u vanishes. Furthermore, the boundary term in the 
associated Brown- York Hamiltonian variation, d(iN(Auj a p A dr/a 13 ) — i^uj a pdrja 13 ), does 
not have a covariant response unless u vanishes. (In that case the reference configuration 
is flat Minkowski space.) From the Komar form of the difference term, we see that their 
quasilocal quantities will have the same values as ours as long as the shift is constant in 
time and the lapse is spatially constant. 

Although our expressions are quite similar to these well known expressions of Brown 
and York, there are some important differences. In particular (i) our expressions are 
manifestly covariant, (ii) we do not require the boundary to be orthogonal to the 
spatial hypersurface (a simplifying restriction they had used which has recently been 
relaxed f)(| |67|, |68|), (iii) we consider more general reference configurations, (iv) our 
displacement vector field N is more general (as will be explained shortly), (v) our 
expressions have a "Komar like" DN term. 



6.1. Static spherically symmetric solutions 

As an example and in order to briefly address the issues of the selection of a reference 
configuration and displacement vector field, we consider certain simple solutions of 
general relativity: the static spherically symmetric metrics. We use the spherical 
coframe 

& = $dt, $ r = $ _1 dr, $ e = rd9, tiP = r smOdcp, (72) 

where $ = $(r). Specifically we have in mind the Schwarzschild (anti)-de Sitter and 
Reissner-Nordstrom type metrics: 



2 



, 2M , 9 Q 
$ 2 = 1 + Ar 2 + . (73) 

In cases like this, where we have an exact analytic form for the metric which depends 
on a few parameters, there is an obvious simple choice for the reference configuration: 
just allow the parameters to have their trivial values. If A ^ we then have the 
choice of a Minkowski or anti-de Sitter reference space. However it is easy to see (with 
computations like those below) that the Minkowski choice will yield an energy which, 
for A 7^ 0, diverges as r — > oo, so we report in detail only the anti-de Sitter choice: 

tf' = $o<ft, ^ r = $o lrfr > °$ e = rde, h = r sin 6 dtp, (74) 
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where $o = 1 + Ar 2 . For the special case of A = the reference space reduces to 
Minkowski space. 

Because of spherical symmetry, linear and angular momentum vanish, so we need 
only calculate the energy. By symmetry, the timelike displacement vector field should 
be orthogonal to the constant t spacelike hypersurfaces. Thus it should have the form 
N = aet, which still allows for different "definitions" of energy. Our preferred choice is to 
define energy using the reference configuration timelike Killing vector d t = $o e t — & e t, 
i.e., a — <E>. Other "obvious" candidates are N = e t and N — e t which correspond to 
the choices of lapse a — 1 and a = $ ~ 1 $, respectively. We will compare the results of 
these choices. 

Under the conditions we are considering, Ar) tr vanishes, in^cJ 3 = OLr\ tc P and 
iNVoP = a&^&oVtoP ■ Consequently 

MN) = ^-Au a p A t]J = -(Au re A Vtr9 + Au r * A Vtrip ) 

X 

= — ( $0 ~ $ )^A^ = — ($ Q -$)rsmOdOAdv, (75) 
X r x 

which leads to the quasilocal energy 

^(iV) = ar ($ -$). (76) 
At the horizon ($ = 0) this reduces to ra$ , giving 

Etiet) = E 6 (dt) = 0, E#(e t ) = r%. (77) 

At large distances E# approaches (a/$ )(M - Q 2 /2r + $Q 2 M 2 /2r). Hence, for A ^ 0, 
the asymptotic result is 

E#(dt) — M — Q 2 /2r, E#(e t ) = E#{e t ) = 0, (78) 

while for A = 

rfl _ m 2 O 2 4- M 2 

E,{e t ) — M — , E4d t ) = E»(e t ) = M — 2r ' (?9) 

On the other hand, for the other control mode 

B U (N) = ^Au% A Q-^wJ = ^(Acu re A ?,tre + Ac/* A j)^) 

= =^ — — - ■&* A ■&* = =^($ -&)r sin 9 d9 A dip , (80) 
X® r x® 

leads to the quasilocal energy 

4(iV) = «rr 1 $ ($ -$). (81) 
At the horizon this is ra($Q/f), giving 

^(e t ) = oo, £^)=r$ 2 = 2(M-Q 2 /2r), = r$ , (82) 
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while for very large r it approaches (M — Q 2 /2r + $q 2 M 2 /2r), giving, for A 7^ 0, 

E u {dt) — M — Q 2 /2r, E u (et) = E u {e t ) = 0, (83) 
and, for A = 0, 

EU°et) = E w (d t )=M-^-^L } E u) {e t ) = M- Li . (84) 

2r 2r 

Because of the vanishing values at large distances for A 7^ 0, e t and e t seem to be 
unsuitable choices. Then the vanishing of E$(dt) at the horizon suggests that analytic 
matching may not be so physical. 

A reasonable alternate way to determine the values of the reference quantities is by 
embedding a neighborhood of the boundary into a space which has the desired reference 
geometry and then pulling back all the quantities to the dynamic spacetime. In the 
present case we can assume that the reference geometry coframe has the same spherical 
form: 

¥ = %dt\ S r = $ _1 (lr' ) ? = rW, ^ = r' sin6' dip', (85) 

where, in particular, $q = 1 + Ar /2 . We first identify corresponding foliations by spacelike 
hyper surf aces. A good criteria for this identification is to relate points that have the 
same trace for the extrinsic curvature. In this particular case this criteria is satisfied by 
the obvious choice: the one form dt' corresponds to f3dt for some function (3. Within 
each spacelike hypersurface a neighborhood of the (assumed to be spherical) boundary 
S = <9£ is diffeomorphic to a suitable image in the reference geometry. A reasonable 
choice for the embedding is to require that the corresponding 2 spheres have the same 



intrinsic geometry. The uniqueness of such an identification has been discussed in f!4 
(One alternative would be to try to match the extrinsic geometries.) In the present 
case it simply means that r = r', 6 = 9', ip = ip'. To completely fix the embedding 
we need to specify the corresponding timelike unit. An obvious simple choice is to take 
t = t', i.e., (3=1. That just gives the analytic matching case already considered. On 
the other hand, giving due consideration to the operational measurement procedure, 
it seems more physically reasonable as well as more geometric to match corresponding 
units of proper time rather than coordinate time. Thus we can choose to fix our relative 

o 

time coordinates by identifying and i?* on <9£, consequently (3 = $$g . It should be 
noted, however, that this choice is not integrable. This type of geometric matching is 
instantaneous. At a later instant of time the procedure will determine a new reference 
configuration. It cannot be expected that these instantaneous reference configurations 
will mesh to form a single reference geometry. 

Proceeding with our quasilocal calculations as before, again, because of spherical 
symmetry we need only calculate the energy and different choices for the displacement 
vector field give different "definitions" of the energy. Our preferred choice is to use the 



24 



reference configuration timelike Killing vector d t > which corresponds in this case to the 
choice of lapse a = $o- Other obvious candidates are N = dt and N = e t = &t which 
correspond respectively to the choices of lapse a = $ and a = 1. We will now compare 
the results of these choices. Under the conditions we are considering, Ar] tr vanishes and 

ijvty/ = arjtoP = ar}to? = InVoP so 

B U = B*= ^Au a /3 A nJ = -(Au re A r) tre + Au r * A rjtr V ) 

= — ( $0 ~ $ )^ g A^ = — ($ Q -$)rfnn9d9Adtp. (86) 
X r x 

Consequently the quasilocal energy is 

E(N) = Etf(N) = E U (N) =ar($ -$). (87) 

For very large r this approaches a($o 1 (M — Q 2 /2r) + $o 3 M 2 /2r) while at the horizon 
($ = 0) it reduces to ar$o- Hence, for asymptotically Minkowski spaces (A = 0, $o = 1) 
our expressions yield 

M 2 -O 2 M 2 + O 2 

E(d t ,) = E(e t ) = M + M - * , E(d t ) = M 



2r v w 2r 

for very large r and 

E(&) = E(e t ) = r, E(d t ) = 0, (89) 

at the horizon. For asymptotically anti-de Sitter spaces our expressions yield 

E{dtf) = E{d t ) =M-Q-, E(e t ) = M$ 1 -> 0, (90) 
for very large r and 

E(dt,)=r& , E(d t ) = 0, E(e t ) = r$ , (91) 
at the horizon. 

We have proposed a reference configuration Killing field as a good choice for the 
displacement vector field. The above calculations show that for geometric matching 
this choice gives reasonable values in practice. Note that the choice N = dt gives a 
vanishing value for the quasilocal energy within the horizon. For asymptotically flat 
space our spherically symmetric quasilocal energy values are the same as those found 
by Brown and York from their Hamilton- Jacobi approach. More generally, our proposal 
for the evolution vector field relaxes their choice of a = 1. Consequently we get the 
total energy M for asymptotically anti de-Sitter solutions whereas their choice of a = 1 
leads to E — M(l + Ar 2 ) -1 which vanishes asymptotically. For the Hamilton- Jacobi 



approach to anti de-Sitter space see [59 



A direct consequence of our favored choice of the reference configuration timelike 
Killing vector is that, just as for Brown and York, the quasilocal energy for the 
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Schwarzschild solution is 2M at the horizon and M at infinity. Thus between the 
horizon and spatial infinity the quasilocal energy is negative, contrary to a property 
favored by several investigators [I7|, ||, [I0| . However, as G. Hay ward |7(J has observed, 
the quasilocal energy cannot be positive in general — simply because closed spaces must 
have zero total energy so they must have negative regions to balance the positive regions. 

The result of Landau-Lifshitz, Tolman [|7TJ and S. Hayward |13[ for the Reissner- 
Nordstrom metric is E = M — Q 2 /2r. They find that the gravitational field has a 
remarkable difference from the electromagnetic field. The energy of the electromagnetic 
field is shared by the interior as well as the exterior of the horizon of the system, but 
the gravitational field energy is confined to its interior only. However, in our result the 
distribution of the gravitational and electromagnetic energy contributions is similar. 
Both are shared by the interior and exterior. One again this is a consequence of our 
choice for N. For the Q = case, some, e.g. |17| , have argued that the quasilocal 



energy for spherically symmetric solutions should simply be M independent of r. In 
our formulation that can be achieved simply by choosing the lapse to have a suitable 
value. (The special orthonormal frame Hamiltonian [[7^] also naturally suggests a choice 
which "localizes" all of the mass within the horizon.) More generally we can obtain the 
aforementioned value E = M — Q 2 /2r simply by choosing a (rather strange) value for 
the lapse such as a = ($ + $)/2- By a judicious choice of the reference configuration 
as well as the lapse our expressions can yield results matching most of the values (e.g., 
those in |[T3||) found by other investigators for spherically symmetric systems. 



7. Discussion 

The quest for gravitational energy-momentum has led to the quasilocal idea: quasilocal 
quantities, including energy-momentum and angular momentum, are associated with 
a closed 2-surface. In this work we have taken a Hamiltonian approach, defining the 
quasilocal quantities in terms of the value of the Hamiltonian for a finite region. This 
value is completely determined by the Hamiltonian boundary term. 

We have considered a rather general class of geometric gravity theories, namely 
those having an independent metric and connection (Metric- Affine gravity). The field 
equations were presumed to follow from a Lagrangian which depends on the metric, 
torsion and curvature. (However, there is no real barrier to extending the procedure to 
include derivatives of these fields.) In our approach, special restrictions including metric 
compatibility, a symmetric connection or teleparallel geometry are easily introduced to 
lead to the familiar theories of interest to most investigators. 

For this general class of geometric gravity theories, by using a covariant canonical 
Hamiltonian formalism, we presented both a general procedure, based on the boundary 
term in the variation of the Hamiltonian and its symplectic structure, and a specific 
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proposal: that this term should have a well defined covariant symplectic structure. We 
found that, for each dynamic field, there is a choice between two covariant boundary 
terms (essentially Dirichlet or Neumann boundary conditions). We then applied our 
formalism to black hole thermodynamics obtaining thereby an expression for the entropy 
in these general theories. Then we restricted our expressions to general relativity and 
applied them to well known spherically symmetric solutions, comparing the results with 
those found by some other authors. (Remarkably one of our expressions turns out to be 
equivalent to an expression recently developed and applied by Katz and his coworkers.) 
An important omission, in the context of Einstein's theory, is the relationship of our 



formalism to spinor formulations. This topic will be considered in a follow up work |74 
Here we wish to further discuss several issues that have arisen. 

From a physical ("operational") point of view, suppose some observers attempt to 
measure the quasilocal quantities. What kind of procedure could they use? What will 
the result depend on? Ideally it should be purely geometric — depending only on the 
dynamical variables at the boundary, independent of any reference frame or choice of 
coordinates. Careful analysis, however, has led us to the recognition that specifying a 
quasilocal gravitational energy-momentum requires many choices, in particular 



i) the theory: e.g., GR, Brans-Dicke, Einstein-Cartan, PGT, MA (Metric- Affine p 



ii) the representation or dynamic variables: e.g., the metric, orthonormal frame, 
connection, spinors; 

iii) the control mode or boundary conditions: e.g., covariant, Dirichlet/Neumann for 
each dynamic variable; 

iv) the reference configuration: e.g., Minkowski; 

v) a displacement vector field: we propose using a Killing field of the reference 
geometry. 

In this investigation we have determined how the quasilocal expressions should depend 
on these choices. 

Given a specific theory and representation our formalism determines the covariant 
quasilocal expressions uniquely, however we do not yet get a unique value for the 
quasilocal quantities in a finite region. This is because there are two features still 
to be determined. One is the displacement vector N, and the other is the reference 
configuration. 

The most ambiguous part of our program is the reference configuration. (We 
prefer the more general term "reference configuration" as opposed to "reference 
geometry" simply because the same type of analysis applies to other dynamic fields, 
e.g., electromagnetism.) It can not be avoided in our expressions. Operationally it 
can be related to the calibration or "renormalization" of our measurement. If the 
dynamical variables take on the reference values our expressions will give zero values for 
all quasilocal quantities. 
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How to choose an appropriate reference configuration? We are not really satisfied 
with our understanding of this topic but will share our current thoughts here. Since our 
boundary could be anywhere we could simple use a globally fixed background geometry. 
For an asymptotically flat or constant curvature spacetime, we can choose the simplest 
global spacetime (let the source vanish in the original spacetime) which has the same 
asymptotic behavior as the physical spacetime. This is not so ambiguous asymptotically; 
it allows one to determine the total conserved quantities. But a global background 
is hardly necessary, moreover it is operationally impractical. How then should one 
determine the reference configuration for a finite region? 

We can imagine the dynamic fields in the neighborhood of the boundary as deviating 
from some "ground state" configuration. These "reference values" can be viewed as 
defining (at least a portion of) a conceptually distinct geometric reference space. A 
slight shift of view then leads to a simple method (which we have elaborated upon) 
for constructing the reference configuration for the geometric quantities of interest 
here: use an embedding (this depends on some matching criteria such as proper 
time and intrinsic 2-surface isometry) of a certain neighborhood of the topologically 
compact orientable boundary spatial 2-surface (usually with spherical topology) into 
a space having the desired reference geometry (ordinarily a Minkowski geometry, but 
alternatives include Schwarzschild, (anti) de Sitter, or a homogeneous cosmology). Pull 
back the geometry to determine the reference configuration. (At the same time we can 
also pull back a displacement vector.) The important question is whether we can find a 
unique embedding? (And, of course, what does it mean physically?) Brown and York 
have proposed embedding the 2-surface isometrically ||14|. They have referred to some 
uniqueness theorems for the embedding of topologically spherical positive curvature 2- 
metrics. Such an embedding seems like a good idea. However, we do not yet adequately 
understand what this means physically or whether some other idea would be equally 
reasonable. 

Another way to deal with the reference configuration is to replace it with some 
additional covariant field satisfying some propagation equation. To our knowledge 



this approach has so far only been used in the context of spinor expressions [[75], [76 



In this case the propagation equation, in effect, implicitly determines the reference 
configuration. For example, the Dougan and Mason spinor propagation equation |J 
determines a spinor field on a closed 2-surface. There is a good prospect for relating 
this approach and the geometric embedding approach. For example a spinor field 
satisfying a suitable propagation equation could be used to determine an orthonormal 
frame which, in turn, could be used to determine a reference configuration for the 
geometrical quantities. 

Turning to the displacement vector iV (which corresponds to the lapse function 
and the shift vector in the ADM analysis |]77f ). How should we choose the appropriate 
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displacement vector in order to get the "proper" quasilocal quantities? For the total 
conserved quantities, the choice is unique, namely an asymptotic Killing vector. But 
for the quasilocal quantities in a finite region, the corresponding Killing vector may not 
exist. Obviously energy is related to time translation, momentum to space translation 
and angular momentum to rotation but which precise choice of timelike vector field 
gives the energy? For example the Brown and York |14| choice for energy, lapse = 1 
and shift = 0, seems natural. Yet with it, for asymptotically anti-de Sitter solutions, 
the quasilocal energy will not converge to the total energy. The choice of N is probably 
best tied to the choice of reference configuration. We have proposed using a Killing 
field of the reference configuration to define quasilocal energy, momentum and angular 
momentum. As in |L4|| alternate choices for the displacement vector field could be used 
to distinguish a quasilocal quantity from a conserved charge. For example, each Killing 
field of the dynamic spacetime would give a conserved charge. 

The quasilocal expressions that we found are well behaved and are related to 
expressions found by other investigators. They satisfy the usual criteria: their values 
have good correspondence limits to flat space and spherically symmetric values, with 
good weak field and asymptotic limits for both asymptotically flat and constant 
curvature spaces, and in both the Bondi and ADM limit. (However a complete analysis 
of the asymptotics which will give well defined quasilocal quantities and the necessary 
conditions for the Hamiltonian to be a differential generator in the sense of [7^ remains 
to be completed.) Although our expressions mesh with a positive total energy proof for 
asymptotically flat solutions to Einstein's theory, our expressions are not locally positive 
in general. However, we regard this as a virtue rather than a drawback for they thus 
allow for the correct (vanishing) total for closed spaces. In any case, such criteria are 
known to be insufficient. 

We have emphasized an additional principle of the Hamiltonian formalism 
concerning the vanishing of the boundary term in the variation of the Hamiltonian 
for suitable boundary conditions. Specifically we have required that the boundary term 
in the Hamiltonian variation have a well defined covariant symplectic structure. 

We stated that we have determined the expressions for quasilocal quantities uniquely. 
Here the meaning of the word uniquely is based on the requirement of the covariant 
symplectic structure. We could release this rather strong requirement to a weaker 
version which only requires a well defined symplectic structure. For example, for a 
vector field W 11 , we could use Dirichlet conditions on part of the spacetime components 
say W°, W 1 and Neumann conditions on the remaining components. Then many more 
quasilocal expressions are possible. This is just the sort of thing that happens with the 
%n factor in the boundary term in the Hamiltonian variations (|T7|JT9|) . 

From the physical point of view, certain noncovariant choices might be distinguished 
and may sometimes be more practical. One approach is to consider the symplectic 
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structure, taking into account the initial value constraints or the boundary conditions 
needed for a good initial value problem. (Boundary control or fixed on the boundary, 
by the way, does not mean having a constant value, but rather means that the value 
of the variable is a preprogrammed function of time.) In particular, for the Einstein 
theory, Jezierski and Kijowski M have investigated decomposing the variables into the 
true unconstrained degrees of freedom. They considered the form of the initial value 
constraints and found that for certain control modes the boundary value problem was 
not well posed. From a naive view point, the symplectic structure reflects the control- 
response relation of the physical system. So for some specific theory, the symplectic 
structure may need some modifications, because the constraints could forbid some 
control mode at the boundary. In principle, we could handle such constraints by 
imposing them with Lagrange multiplier terms in the variational principle. This may 
modify the symplectic structure and the expressions for quasilocal quantities, effectively 
automatically including the constraints. Further investigation is certainly needed on this 
topic. 

Others boundary control choices might be more closely related to what an observer 
would directly measure. For example, in an electromagnetic system, we can naturally 
fix the electric potential and the tangential components of the magnetic field instead 
of the whole 4-vector potential at the boundary, or we can fix the normal component 
of the electric field and the tangential components of the vector potential (essentially 
the normal component of the magnetic field). The physical meaning as well as the 
procedure for implementing these choices is well understood. But the boundary term in 
the corresponding Hamiltonian variation for such control modes depends on the fields in 
a noncovariant way. However, these non-covariant expressions are related to a covariant 
expression of our type by a simple Legendre transformation on the boundary. Having 
measured the values for one of these situations, we could calculate the value of our 
covariant quasilocal expressions even if they could not be conveniently measured directly. 

We have emphasized the importance of covariance for good physical quasilocal 
quantities. Yet here exists a noncovariant term, i^u^ApJ 3 , in our Hamiltonian 
boundary expressions. This term merits deeper consideration. We noted that it was 
through this term that our expressions mainly differed from those found by others. 
Consequently, further investigation is desirable to ascertain its importance and how 
essential it is that it have exactly the form that we found. We have decomposed this term 
by using the identity (24). The first two terms on the rhs of this identity are covariant 
and represent a real physical effect. The last term, which reflects a noncovariant property 
of our expressions, is an unphysical term presenting a dynamical reference frame effect. 
This situation is like the centrifugal force in mechanics. The complete Hamiltonian 
generates the (arbitrary) dynamical evolution of the frame (effectively, the observer) 
as well as the physical variables. Through the relation (E4T), the complete Hamiltonian 



30 



implicitly includes a time derivative term on the boundary. Within the context of 
Einstein's theory, boundary time derivative terms have been noted some time ago by 
Kijowski [^7| and more recently by G. Hayward ||7D| . 

A major virtue of our formalism is that it provides for system and order. Although 
it allows for many "energies", yet each has (in principle) a clear physical meaning. 
The formalism systematically associates a quasilocal energy-momentum with a specific 
Hamiltonian boundary term (determined by the choice of theory, representation — 
dynamic variables, control mode, reference configuration, matching criteria, and 
displacement vector field). The gravitational quasilocal energy-momentum is thus, 
like the energies in thermodynamics, connected with a definite physical situation 
with definite conditions on the boundary. But much more investigation is needed to 
understand the significance of various boundary conditions. Hence our formalism shifts 
the attention to questions like "What is the real physical significance of controlling the 
metric (or the connection) on the boundary of a finite region? Which boundary control 
(if any) is the analog of thermally insulating a system? Which corresponds to a thermal 
bath?" Some insight is afforded by classical electrodynamics where the answers to such 
questions are understood |27|, ||. 

Of course one way to further investigate the various quasilocal expressions is to do 



more direct calculations for exact solutions, e.g., fT5| . However, a deeper theoretical 
investigation should be more revealing. Our formulation provides a good starting 
point for such an investigation. Note that all of the expressions presented here 
correspond to the work done in some (ideal) physical process. The situation is similar 
to thermodynamics with its different energies (enthalpy, Gibbs, Helmholtz, etc.) An 
even better analogy is the electrostatic work required while controlling the potential on 
the boundary of a region vs. that required while controlling the charge density . 

However, we have gone far enough for the present. Our intention here was primarily 
to set up a framework, to indicate something of its scope and the efficacy of the 
principles. In follow up works we plan to extend these investigations in several directions 
including the application to and relationship with spinor expressions for quasilocal 
quantities |74] . 
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